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3 : ABSTRACT 

o ■ 

^ , We examine the coarse-grained phase-space density profiles of a set of recent, hi^- 

, ^ \ resolution simulations of galaxy-sized Cold Dark Matter (CDM) halos. Over two and a 

^ I half decades in radius the phase-space density closely follows a power-law, p/a"^ cx r~", 

I ' with a ~ 1.875. This behaviour closely matches the self-similar solution obtained by 

■ Bertschinger for secondary infall of gas onto a point-mass perturber in a uniformly ex- 
T-H ■ panding universe. On the other hand, the density profile corresponding to Bertschinger's 

, solution (a power-law of slope r'^°^~^) differs significantly from the density profiles of 

i CDM halos. CDM halo density profiles are clearly not power laws, and have logarith- 

\ mic slopes that gradually steepen with radius, roughly as described by Navarro, Frenk 

' &: White (NEW). We show that isotropic, spherically-symmetric equilibrium mass dis- 

' tributions with power-law phase-space density profiles form a one-parameter family of 

' structures controlled by the ratio of the local velocity dispersion to the "natural" veloc- 

Q_i' ity dispersion at some fiducial radius, ro; k = 47rG/>(ro)rQ/(T(ro)^. Eor k = a = 1.875 

Q . one recovers the power-law solution p cx r^"~^. As k increases, the density profiles 

^ , become quite complex but still diverge like r'^""^ near the center. Eor k larger than 

\ some critical value, Kcrit(a)i solutions become non-physical, leading to negative den- 

^ \ sities near the center. The critical solution, k = Kcriti corresponds to the case where 

' the phase-space density distribution is the narrowest compatible with the power-law 

• phase-space density stratification constraint. Over three decades in radius the criti- 

■ cal solution is indistinguishable from an NEW profile, although its logarithmic slope 
asymptotically approaches — 2a/5 = —0.75 (rather than —1) at very small radii. Our 
results thus suggest that the NEW profile is the result of a hierarchical assembly process 
that preserves the phase-space stratification of Bertschinger's spherical infall model but 
which "mixes" the system maximally, perhaps as a result of repeated merging, leading 
to a relatively uniform phase-space density distribution across the system. This finding 
offers intriguing clues as to the origin of the similarity in the structure of dark matter 
halos formed in hierarchically clustering universes. 

Subject headings: cosmology: dark matter — cosmology: theory — galaxies: formation 
— galaxies: structure — methods: analytical — methods: numerical 



^CIAR Scholar and Alfred P. Sloan Research Fellow 



- 2 - 



1. Introduction 

In the currently favored paradigm for the formation of structure in the Universe, luminous 
objects such as galaxies and galaxy clusters are embedded within extended hales of cold dark matter 
(CDM). The non- linear equilibrium structure of these halos has long been thought to contain clues 
as to the astrophysical nature of dark matter, a fact that has fueled strong interest in detailed 
theoretical predictions for the structure of dark matter halos, as well as in constraints placed by 
observations of disk galaxy rotation curves (Frenk et al. 1988; Flores et al. 1993; Flores & Primack 
1994; Moore 1994; Moore et al. 1999b, McGaugh k de Blok 1998, van den Bosch et al. 2000; Swaters, 
Madore & Trewhella 2000; van den Bosch & Swaters 2000), by gravitational lensing of galaxies and 
clusters (see, e.g., Tyson, Kochanski & dell' Antonio 1998, Williams, Navarro & Bartelmann 1999), 
and by detailed studies of the stellar dynamics of spheroidal galaxies (Carollo et al. 1995; Rix et 
al. 1997; Gerhard et al. 1998; Cretton et al. 2000; Kronawitter et al. 2000). 

Unfortunately, observational constraints tend to be strongest near the center of dark mat- 
ter halos - where most of the luminous material in galaxies resides - but also where theoretical 
predictions are least robust. This is because of the difficulties inherent to simulating accurately 
through N-body methods regions where overdensities exceed ~ 10*^ and where particles may have 
completed thousands of orbits during a Hubble time. In spite of these difficulties, there has been 
steady progress in our theoretical understanding of the non-linear structure of virialized dark mat- 
ter halos, driven largely by direct numerical simulation. Following on pioneering numerical studies 
by Quinn, Salmon & Zurek (1986), Frenk et al. (1988), Dubinski &: Carlberg (1991), and Crone, 
Evrard k Richstone (1994), Navarro, Frenk & White (1996, 1997, hereafter NFW) found that, 
regardless of mass and of the value of the cosmological parameters, the density profiles of dark 
matter halos formed in various hierarchical clustering cosmogonies were strikingly similar. This 
"universal" structure can be characterized by a spherically-averaged density profile which differs 
substantially from the simple power-laws, p{r) oc r~^, predicted by early theoretical studies (Gunn 
& Gott 1972, Fillmore & Goldreich 1984, Hoffmann & Shaham 1985, White & Zaritsky 1992). The 
profile steepens monotonically with radius, with logarithmic slopes shallower than isothermal (i.e. 
P < 2) near the center, but steeper than isothermal {f5 > 2) near the virial radius of the system. 

The NFW result has been confirmed by a number of subsequent studies (Cole k Lacey 1996, 
Huss, Jain k Steinmetz 1999, Fukushige k Makino 1997, Moore et al. 1998, Jing k Suto 2000), 
although there is some disagreement amongst authors regarding the innermost value of the loga- 
rithmic slope. NFW argued that a fitting formula where l3 = {1 + 3y)/(l -|- y) (where y = r/vg is 
the radius in units of a suitably defined "scale radius" r^) provides a very good fit to the density 
profiles of simulated halos over two decades in radius. Moore et al. (1998), Ghigna et al. (2000), 
and Fukushige k Makino (2000) have argued that (3 converges to a value of ~ 1.5 near the center, 
rather than the 1 expected from the NFW fit. Kravtsov et al. (1998) initially obtained much shal- 
lower inner slopes (/? ~ 0.7) in their numerical simulations, but have now revised their conclusions; 
these authors now argue that CDM halos have steeply divergent density profiles but, depending 
on evolutionary details, the slope of a galaxy-sized halo at the innermost resolved radius may very 
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between -1.0 and -1.5 (Klypin et al. 2000). 

Steep inner slopes have been traditionally disfavored by rotation curve data, a fact that has 
often been used to "rule out" CDM as a viable cosmogony (Flores k, Primack 1994, Moore 1994, 
McGaugh & de Blok 1998, Blais-Ouellette et al. 2001, Cote ct al. 2000, de Blok et al. 2001). 
However, recent reanalysis of the available data suggests that most rotation curves are broadly 
consistent with "cuspy" dark matter cores, provided that the innermost slope is shallower than 
(3 ~ 1.5 and that the "concentration" of the halos is low (Navarro 1998, van den Bosch et al. 2000, 
van den Bosch k, Swaters 2000). 

Although there may not be at present broad consensus regarding how steep the innermost slope 
is, or even whether there should be a well-defined asymptotic innermost slope, there is agreement 
that it will take extraordinary computational effort to reach a robust resolution of the controversy. 
What is required is a statistically significant sample of galaxy-sized halos with sub-kpc resolution, 
an extremely onerous computational task that will stretch the capabilities of the most powerful 
massively parallel computers. Steps in this direction are currently being taken (e.g.. Power et al., 
in preparation), but it will take some time until these efforts yield conclusive results. 

From the theoretical point of view, a number of plausible argTimcnts have been advanced 
in order to try and explain the innermost behaviour of dark matter density profiles from stellar 
dynamical principles. These efforts, however, tend to give non-unique results and have so far been 
unable to explain the remarkable similarity in the structure of dark matter halos of widely different 
mass formed in a variety of cosmogonies (Evans &; Collett 1997, Syer & White 1998, Nusser & Sheth 
1999, Lokas & Hoffman 2000). In this paper we investigate an empirical alternative to analytic 
efforts addressed at estimating the innermost slope of the density profile. Our proposal exploits 
the similarity between the phase-space density profiles of CDM halos and that of the self-similar 
solution for spherical collapse in an expanding universe found by Bertschinger (1985). In addition, 
this offers an attractive scenario for understanding the shape of halo density profiles as well as a 
powerful tool for estimating their slope near the center. 

The outline of this paper is as follows. In Section 2, we present the phase-space density profiles 
for galaxy-sized CDM halos in several recent simulations, and show that they are well approximated, 
for over two decades in radius, by a power law. In Section 3, we investigate equilibrium density 
profiles consistent with this constraint by solving the Jeans equation for spherically symmetric 
systems with isotropic velocity dispersion tensors. Section 4 compares these results with the results 
of numerical simulations. We discuss and summarize our findings in Section 5. 

2. The Phase-Space Density Profile of CDM Halos 

We have chosen for our analysis three high-resolution simulations of galaxy-sized CDM halos 
identified at z ~ 0. Each of these simulations has of order 10^ particles within the virial radius of the 
halo, and was run with PKDGRAV, a massively-parallel N-body code developed by Joachim Stadel 
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and Thomas Quinn at the University of Washington. These simulations are amongst the largest 
carried out so far for galaxy-sized objects, and represent a major investment in computational 
resources. 

Two of the halos were run in the former "standard" CDM cosmology (SCDM, = 1,A = 
0,h = 0.5, as = 0.7) and have circular velocities of ~ 180 and 160 km s~^, respectively. These 
circular velocities are measured at the "virial" radius, r2oo, where the mean inner density of the 
system is 200 times the critical value for closure. These two halos are part of the "Local Group" 
simulation reported in Moore et al. (1999a). The third halo has a circular velocity of 200 km 
and was run in the currently popular LCDM cosmogony (Q^ = 0.3, A = 0.7,^ = 0.65, ag = 0.9). 
Gravitational softenings were chosen to be ~ 3 x 10~^ and 1.5 x 10^"^ r2005 for the SCDM and 
LCDM runs, respectively. PKDGRAV uses a multi-stepping algorithm to integrate the equations 
of motion; particles with the smallest time bins may take up to 100, 000 timesteps to evolve from 
the initial redshift {zi = 99 for SCDM and Zi = 49 for LCDM) to the present. The initial conditions 
for the SCDM runs were setup using the algorithms devised by the N-body shop at the University 
of Washington. Those for the LCDM run were set up in a completely independent way using the 
algorithms described by Navarro, Frenk & White (1996, 1997; see also Efstathiou et al. 1985). As 
we describe below, the results of these two independent runs are consistent with each other, which 
suggests that our conclusions are independent of both the cosmological model adopted and of the 
particular choice of initial conditions setup algorithms. 

Figure 1 shows the spherically-averaged phase-space density profile of the three CDM halos. 
Solid (dashed) lines are used for the SCDM (LCDM) halos. The phase-space density is computed 
in spherical bins containing 2, 000 particles each and is defined as p/a^, where p is the mass density 
and a is the 1-D velocity dispersion in the bin.^ For ease of comparison, we have chosen to 
normalize p/cr'^ so that all three curves coincide at 0.01 r2oo- The important point illustrated by 
this figure is that, over more than two decades in radius, the phase-space density profile is very well 
approximated by a power law of slope —1.875 (thin straight line). This is quite remarkable, given 
that both the density profiles (shown in Figure 2) and the velocity dispersion profiles of these halos 
deviate quite strongly from simple power laws, as described by NFW. 

Also remarkable is that the slope of this power-law coincides with the self-similar solution 
derived by Bertschinger (1985) for secondary infall onto a spherical perturbation in an unperturbed 
Einstein-de Sitter universe: Bertschinger 's solution is plotted with solid circles in Figure 1. This 
solution corresponds to the self-similar equilibrium configuration of a 7 = 5/3 gas formed by 
spherical accretion onto a point-mass perturber in an otherwise uniform Einstein-de Sitter universe. 
The quantity shown by the solid circles is the quantity "equivalent" to the phase-space density, 
p5/2^p3/2^ where P is the local (isotropic) pressure. As discussed by Bertschinger, this solution is 
the most appropriate to compare with our numerical results for CDM halos, given that the velocity 



^The velocity dispersion tensor is roughly isotropic near the center and only mildly radially biased in the outer 
regions. 
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dispersion tensor in this case is only mildly anisotropic. Radii are normalized assuming that r2oo 
equals the shock radius of the self-similar solution, which corresponds to roughly one-third of the 
turnaround radius. The vertical normalization is arbitrary and has been chosen to match the N- 
body results at 0.01 r2oo- Taking p5/2/p3/2 ^ (T/p"'^^)^^/'^"' to be a measure of the local "entropy" 
of the system, Figure 1 shows that CDM halos have the same radial entropy stratification as the 
simple spherical collapse solution. It is possible that this power-law stratification is a fundamental 
property which underlies the similarity of structure of cold dark matter halos. 



3. Density Profiles 

Density profiles consistent with the power-law phase-space density profile shown in Figure 
1 can be obtained by assuming hydrostatic equilibrium. For an isotropic, spherically-symmetric 
system of coUisionless particles, the Jeans equation may be written as 

d{pa^__ d^__ GM{< r) 
dr ar 

where $ is the gravitational potential, and M{< r) is the mass interior to r (Binney and Tremaine, 
1987, p. 198). Equation (1) is equivalent to the equation of hydrostatic equilibrium for a gas of 
pressure P = pa^. Dividing both sides of the equation by —Gp/r, and taking derivatives with 
respect to r, we can rewrite this equation as, 

' / ^'^ ^ M _ .M{< r) = Anpr^, (2) 



dr \ Gp \ dr J J dr 

where the last equivalence applies to a self-gravitating system. Assuming that the phase-space 
density is a power-law of radius, i.e. that 

p/a\r) = ipo/4)ir/ro)-", (3) 

where tq is an (arbitrary) reference radius and ctq = o'(^o) ^-nd po = p{ro), and defining the 
dimensionless variables, x = r/ro and y = p/po, (2) can then be written as, 

where 

K^l^ (5) 

is a dimensionless measure of the velocity dispersion at tq. A simple interpretation of the parameter 
K can be found by choosing ro to be rp, the radius where the circular velocity of the system peaks. 
In this case, and since p{rp) = M(< rp)/4'Krp, we can rewrite k as, 

_ AnGpprl _ GM{< rp)/rp _ 

<jp yjp <Jp 
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In other words, k measures the local velocity dispersion of the system in units of its circular velocity 
at the radius where the circular velocity peaks. 

Equation (4) has the form, 

y" + f{y',y,x) = 0, (7) 

(where y' = dy/dx,y" = (fy/dx^) so given initial conditions y{l) = p{ro)/po = 1 and y'{l) = y'l 
at X = l(r = ro), we can solve it numerically for r > ro or r < tq by integrating inwards or 
outwards. The resulting solutions are parameterized by a, y[ and k, so for a = 1.875, wc obtain a 
two-parameter family of solutions. Clearly, integrating (7) with respect to some different variable 
z = tx, where t is a constant, will produce an identical solution, provided that we scale y[ by t and 
K by t2a/3-2^ gQ £g^g^ ^Yie solutions are degenerate and there is effectively a single free parameter 
(which we can take to be k) that characterizes the full solution set.^ 

Equation (7) admits power-law solutions. Solving for y = x~^ gives the following conditions. 
For K = we have, 

(5 = 2a/5 

and 13 = a -h 3/2, (8) 

while for K 7^ the conditions are, 

13 = 6 -2a 

K = 8(a -3/2)(5/2 - a). (9) 

For example, the singular isothermal sphere is recovered for k = a = /? = 2. 

Thus assuming k ^ 0, that is, a finite velocity dispersion at ro, there is a single power- 
law solution to (7) with slope /3 = 6 — 2a = 2.25 for a = k = 1.875.^ It is convenient to set 
= 2a — 6 = —2.25 and to generate the family of solutions simply by varying k. Figure 2 shows 
a few of these solutions (for a = 1.875), labeled by the particular value of k beside each curve. 
With the exception of the power law, the shape of the equilibrium density profiles is complex, with 
an outer cutoff, one or more inflection points, and an inner cusp. For n larger than some critical 
value, /Ccrit(a), the solutions become non-monotonic, and the density vanishes at a finite radius. 
All solutions with 1.875 < k < Kcrit appear to have a steep inner cusp with asymptotic slope close 
to /? = 6 — 2a, but the critical solution approaches asymptotically (3 = 2a/5 = 0.75 as r tends to 0. 



^This assumption breaks down if we choose y[ to be positive, or negative but too small. In this case p(r) increases 
with r and the solutions are unstable and lack direct physical interpretation. 

■* Curiously, besides the singular isothermal sphere this is the only other case where the power-law solution has 
K = a. 
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4. Comparison with Numerical Density Profiles 

As discussed in the previous section, the isotropic Jeans equation admits a family of solutions 
for the density profile under the constraint pja^ oc r~". The family includes a power-law, p cx r~^, 
with /3 = 6 — 2a = 9/4 = 2.25 (for a = 15/8 = 1.875), which corresponds to Bertschinger's spherical 
infall solution. This is shown with filled circles in Figure 3, where wc also show the density profiles 
corresponding to the CDM halos in the N-body simulations. Radii have been scaled here, as in 
Figure 2, to the radius where the slope of the density profile is —2.25, and densities to the critical 
density. Clearly, the power-law solution is a poor fit to the result of the N-body simulations, which 
become noticeably shallower than r""^/^ near the center. 

The N-body results can be well fitted by the profile proposed by NFW, which is shown by 
the open triangles in Figure 3. Interestingly, over approximately three decades in radius, the shape 
of the NFW profile is essentially indistinguishable from the "critical" solution alluded to above. 
Discrepancies occur only at large radii, where substructure leads to systematic deviations from 
the power-law behaviour for the phase-space density. The meaning of the critical solution may 
be illustrated by considering the phase-space density distribution corresponding to the different 
solutions. This is shown in Figure 4 for the solutions illustrated in Figure 2, for systems normalized 
to have the same total mass, energy, and mean density inside some fiducial radius. Clearly, the 
phase-space distribution function is broadest for the power-law solution {k = a = 1.875) and gets 
increasingly narrower as k increases; the critical solution (k = Kent) corresponds then to the most 
sharply peaked phase space distribution compatible with a monotonically decreasing (non-hollow) 
density profile and with the entropy stratification constraint. The critical solution may thus be 
interpreted as a "maximally mixed" configuration where the phase-space density is as uniform as 
possible across the system. 

This leads to the following interpretation for the origin of the NFW profile. Gravitational 
assembly of CDM halos leads to a simple power-law radial stratification of the phase-space density. 
If spherical symmetry is imposed, as in the case treated by Bertschinger (1985), the collapsing 
radial mass shells generate progressively larger "entropies" (lower phase-space densities) as they 
pass through the shock and settle into hydrostatic equilibrium, leading to steeply cusped power-law 
profiles with slope f3 = 6 — 2a. On the other hand, when the assumption of spherical symmetry 
is released and the collapse proceeds through many stages of hierarchical merging, mass shells are 
continuously "mixed" and the profiles tend to the critical solution: that corresponding to the most 
uniform entropy distribution compatible with a non-hollow density profile and with the power-law 
entropy stratification constraint. This density profile closely resembles the NFW profile over a large 
dynamic range in radius. 

Finally, it is important to note that, despite the similarity shown in Figure 3, there are impor- 
tant differences between the NFW profile and the critical solution, the most notable being that the 
latter tends to an asymptotic central slope of —f3 = —2a/ 5 = —0.75 rather than to —1 as in NFW's 
fitting formula. A simple approximation to the radial dependence of the slope of the "critical" 
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density profile is given by 



dlnx 



dlnp 



0.75 + 2.625 x^/^ 
l + 0.5xV2 



(10) 



which is accurate to 3% for x = r/ro < 4. Here tq is the radius where the logarithmic slope of the 
density profile equals —2.25 and ro = (5/3) r^, in terms of the NFW scale radius Vg- 



The power-law phase-space density profile thus offers a natural way to describe the structure 
of dark matter halos at radii where simulations become increasingly difficult and expensive, but 
also where observational constraints are strongest. Provided that the velocity dispersion tensor 
remains nearly isotropic, the "critical" solution provides a clear prediction as to the behaviour 
of the logarithmic slope of the density profile: it should become progressively shallower towards 
the center, converging asymptotically to a value of —2a/5 = —0.75. This is interesting since, as 
mentioned in §1, slopes shallower than —1.5 appear to be consistent with the recent reanalysis of 
the rotation curve dataset by van den Bosch et al. (2000) and by van den Bosch & Swaters (2000). 
It is important to stress, however, that a shallow central slope does not guarantee consistency 
with observations, which constrain the detailed radial dependence of the density profile slope as 
well. Our results do, however, offer a clear prediction for extrapolating the mass profiles to regions 
that are very difficult to probe numerically. Is the extrapolation of the power-law behaviour to 
very small radii warranted? This question ultimately will have to be answered by direct numerical 
simulation, although there is no obvious a priori reason why a power law scaling that is valid for 
over two decades in radii should break down nearer the center. 

One important point to note is that the critical solution is clearly at odds with the proposal of 
Moore et al. (1998) and Ghigna et al. (2000) that the innermost slopes of CDM halos converge to a 
value not shallower than about —1.5. However, it should be emphasized that their conclusion was 
based on the simulation of a single halo simulated in a standard CDM universe (SCDM) and on a 
significantly different mass regime (galaxy clusters) than probed here. Thus the possibility remains 
that this particular system may not be representative of the general population or that the density 
profiles of clusters are steeper than those of galaxy-sized halos. Since these authors use the same N- 
body code as in the present work and, indeed, we use their own results for galaxy-sized halos in this 
manuscript, it is unlikely that the discrepancy is due to subtle errors associated with the numerical 
setup of the simulations. Confirming which of these possibilities holds will require a statistically 
significant sample of halos simulated with resolution comparable to the systems used here. Finally, 
it is also possible that the cluster simulated by Moore et al. (1998) and Ghigna et al. (2000) differs 
from the galaxy-sized halos we present here in other, more subtle ways. For example, it may be 
significantly more triaxial than the systems analyzed here, or perhaps its velocity dispersion tensor 
is very anisotropic, in conflict with the assumptions of this work. Again, a detailed reanalysis of 
the discrepant system, extended to a statistically meaningful sample, appears necessary in order 



5. 



Discussion 
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to explain this discrepancy conclusively. We are currently working on this issue and plan to report 
our findings soon (Stadel et al., in preparation). 

6. Summary 

We examine the spherically averaged phase-space density profiles of N-body simulations of 
CDM halos in the SCDM and LCDM cosmologies and find that they are very well approximated 
by a power-law {p/a^ oc r"^'^''^) over more than two decades in radius. The slope of this power 
law is consistent with that of the self-similar solution for spherical secondary infall derived by 
Bertschinger (1985). Assuming that the phase-space density profile is a power law, and assuming 
isotropy, hydrostatic equilibrium is satisfied by a family of density profiles controlled by a single- 
parameter. This parameter is determined by the ratio of the velocity dispersion to the "natural" 
velocity dispersion of the system at some fiducial radius, tq, k = 47rGp(ro)rQ/cr(ro)^. The parameter 
K can also be expressed in terms of the ratio of the circular velocity to the 1-D velocity dispersion 
at Tp, the radius where the circular velocity peaks {n = F(rp)^/c7(rp)^). For k = 1.875, the density 
profile is a power-law, which agrees with Bertschinger's solution. As k increases the density profiles 
become increasingly curved, although they still approach the steep power-law divergent behaviour 
near the center. For k greater than some criticcil value Kcrit 

~ 2.678 the density profiles become 

unphysical, vanishing at some finite radius near the center. 

The "critical" density profile (k = Kait) corresponds to the maximum value of k consistent 
with a non-vanishing density profile at the center and is essentially indistinguishable from the pro- 
file proposed by Navarro, Frenk & White (1996, 1997) for over two decades in radius. This solution 
corresponds to the most uniform (more sharply peaked) phase-space distribution compatible with 
a monotonically decreasing density profile and with the entropy stratification constraint. This 
suggests that the structure of CDM halos is determined by a radial phase-space density stratifi- 
cation process similar to that established through collapse onto a point mass perturbation in an 
unperturbed expanding universe, and by the uniformization of phase-space density that occurs 
presumably as a result of the many merger and satellite accretion events that characterize the 
assembly of a CDM halo. This identification leaves a couple of important questions unanswered, 
however: (i) why should the phase-space density be a power law of radius?, and (ii) why is the 
exponent the same as in Bertschinger's self-similar solution?. Although we have no clear answer 
to these questions at this point, our results suggest that explaining the origin of the structural 
similarity of CDM halos pointed out by Navarro, Frenk &: White may entail unraveling why the 
radial stratification of phase-space density in CDM halos is the same power-law generated by the 
simple spherical collapse model. Discovering a mechanism that achieves this may provide a simple 
explanation for the universal structure of cold dark matter halos. 

The authors wish to thank Joachim Stadel and Tom Quinn for kindly making their massively 
parallel code PKDGRAV available as well as for their help and guidance running the simulations. 
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Fig. 1. — The phase-space density profiles of three galaxy-sized CDM halos. Solid lines correspond 
to the SCDM halos and dashed lines correspond to the LCDM halo. Vertical normalizations are 
arbitrary and have been chosen so that the curves coincide at about 0.01 ?'2oo- Radii are normalized 
to the virial radius, ^200- The solid circles indicate the self-similar solution obtained by Bertschinger 
(1985) for spherical infall of gas onto a point mass perturber in a uniform Einstein-de Sitter universe. 
Radii for this solution have been normalized by assuming that the shock radius in the solution equals 
r2oo- A power-law of slope —1.875 is shown for comparison (thin solid line). 
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Fig. 2. — Density profile solutions to the isotropic, hydrostatic equilibrium (Jeans) equation, for 
spherically symmetric systems with a power-law phase-space density, p/cr^ oc r~^'^'^^. The curves 
are generated by setting /3 = — d(logp)/d(logr) = 2.25 and p(ro) = po = 1 at some fiducial radius, 
ro, and by varying the free parameter k. The value of k is used to label each curve. 
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Fig. 3. — Density profiles of the CDM halos and the Bertschinger solution, compared with the 
analytic solutions shown in Figure 2. These profiles deviate significantly from the power-law solution 
(p oc r~^'^^) and from Bertschinger's solutions, but are well fitted by an NFW profile (open triangles) 
and by the "critical" (k = Kcrit) solution (see text). 
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Fig. 4. — The phase-space distribution function of the five solutions shown in Figure 2, normahzed 
so that ah have the same total mass, energy, and mean density. The straight line corresponds 
to the power-law solution, k = 1.875. As k increases the distribution becomes more sharply 
peaked. Dashed lines correspond to the hollow (non-physical) density profiles. The critical solution, 
1^ = ^crit = 2.678, corresponds to the most sharply peaked phase-space density distribution or, 
equivalently, to the case where phase-space density is most uniformly distributed across the system. 



